Majorana bound state in rotating superfluid ^He-A between parallel plates 
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A concrete and experimentally feasible example for testing the putative Majorana zero energy state 
bound in a vortex is theoretically proposed for a parallel plate geometry of superfluid ^He-A phase. 
We examine the experimental setup in connection with ongoing rotating cryostat experiments. The 
theoretical analysis is based on the well-established Ginzburg-Landau functional, supplemented by 
microscopic calculations of the Bogoliubov-de Gennes equation, both of which allow the precise 
location of the parameter regions of the Majorana state to be found in realistic situations. 
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Much attention has been focused on Majorana zero- 
energy states in various research fields[H, for example, 
leptogenesis in cosmology in connection with the matter- 
antimatter imbalance problem the fractional quantum 
Hall state, chiral superconductors, and p-wave superflu- 
ids in neutral atomic gases. Of particular interest is the 
possible application to quantum computing, based on 
the fact that a pair of Majorana states is intrinsically 
entangled and topologically protected from external dis- 
turbance. For chiral superconductors and p-wave super- 
fluids, the key is finding the Majorana zero energy state 
bound in a vortex. Its existence is guaranteed by the in- 
dex theorem[lj, implying in this context that in spinless 
chiral p-wave pairing p^ ± ipy of two dimension (2D), an 
odd winding number vortex always possesses a Majorana 
state at exactly zero energy relative to the Fermi level. 
This state is localized at the coreQ. 

To date there is has been no experimental evidence 
to prove its existence in either chiral superconductors Q 
or p-wave superfluids in neutral atomic g neither 
of these phenomena have yet been realized. Therefore, 
concrete examples of p-wave superfluidity are needed. In 
this respect, the ■^He-A phase is a prime candidate for 
testing the Majorana state. 

The ABM (A) phase in superfluid ^He is character- 
ized by the order parameter (OP) = Aoc?^(n -|- 
irn)i (i — x,y,z). The spin and orbital structure 
of p-wave pairing, which is generally defined as A = 
* S/ii (A^ia^pi) (Jy, are characterized by the d vector and 
I vector, respectively!^ [^. (l = ft x rn, rh and ft are unit 
vectors, forming a triad.) It is known that the I vector al- 
ways points perpendicular to a wall surface and d vector 
is perpendicular to the applied field direction if its mag- 
nitude is sufficiently strong {H > Hd — 2.0 mT, where 
Hd is the dipole field[^). Therefore, we can control the 
pairing symmetry both by the confining geometry and 
the field. This kind of controllability in superfluid '^He is 
a great advantage over other research fields, in addition 
to the decisive fact that the precise OP forms of the p- 
wave pairing are established for the ABM (A) and BW 
(B) phasesdS]. 



We consider a parallel plate geometry with a gap com- 
parable to the dipole coherent length ~ 10 /im and a 
field [H > Hd) applied along the plates. In this configu- 
ration the two vectors I \\ d are locked together perpen- 
dicular to the plates to maximize the dipole and mag- 
netic energies, and hence the OP is A^i = Aod^{n + irn)i 
{i — x,y) with /x being fixed. This is a "spinless" chiral 
superfluid, a situation ideal for the Majorana zero energy 
state[l[. 

This parallel plate geometry has been already realized 
0, S|, namely a '^He-A phase sample'^ confined in a 
thin cylindrical region with thickness 12.5 /um and ra- 
dius R — 1.5 mm. This thin plate with a cylindrical 
region is stacked to form a layered structure consisting 
of 110 plates. In an effort to observe a half-quantum 
vortex, the structure was rotated up to a rotation speed 
of 51=6.28 radian/sec (rad/sec) in a rotating cryostat at 
ISSP, University of Tokyo. The texture transition under 
rotation was monitored by an NMR spectrum to identify 
each vortex. 

The experimental geometry of this previous experi- 
ment is ideal for realizing the Majorana zero energy state, 
and hence we have adopted it for our study. The neces- 
sary conditions for the realization of the Majorana state 
are: (1) The bulk superfluid has a chiral p-wave symme- 
try. (2) The vortex is singular with an odd winding num- 
ber. Condition (1) is satisfied because the '^He-A phase 
is orbitally triplet. As mentioned, the d and I vectors are 
locked together perpendicular to the plane. The super- 
fluid state is in the chiral p-wave symmetry characterized 
by P± ~ l^iPx i ipy)/\/2 with frozen spin freedom. 

However, we must carefully examine condition (2). 
There are several possible vortex states under this con- 
figuration, either singular or non-singular vortices, and 
compete with each other energetically. The Majorana 
zero energy state is supported only by the singular vortex 
in general. Thus it is important to determine the possible 
stable parameter space spanned by the rotation speed f2, 
temperature T, radius R of the system, andjDressure P 
in terms of realistic experimental conditions IS]. We also 
study the microscopic excitation spectrum of the Majo- 
rana state on the three-dimensional (3D) Fermi surface 
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by solving the microscopic Bogoliubov-de Gennes (BdG) 
equation for a vortex. Previous theoretical studies have 
been limited to strictly 2D systems [I]- This study clar- 
ifies the detailed spectral structures in the bulk and at 
the boundary where the edge current flows, i.e., the dis- 
persion relation of the Majorana zero mode in 3D. 

We start with the following well-established Ginzburg- 
Landau (GL) functional |5!, |6( to examine the stable re- 
gions for the Majorana state. This GL form is known 
to be adequate for describing the ^He-A phase where the 
OP is described by Azi = dzAi: 

f — fbulk + fgrad (1) 

hulk = -aA*A, + PuA*A*A,A, + f32i5A* A*A,A, (2) 

fgrad = 'K [A*d^d,Aj + A*dAAj + A*d^d,A,] (3) 

where di = Vi — x r)i, is the mass of a '^He 

atom, K is the effective mass, a = aoCL — T/Tc), Pis = 
I3i + and /3245 = /?2 + /34 + /^sd, Hi- We introduce 
A± =t{Ao,T iAy)/V2 with A = A,p, = A+p+ + Aqpq + 
A-P-, {pz = po)- Here we take into account terms up to 
the fourth order invariants and drop the spin degrees of 
freedom. As mentioned, the parallel plate geometry has 
a narrow gap such that the /-vector is fixed normal to the 
plate surfaces, implying that Az = Q- 

The free energy is minimized numerically for var- 
ious radii R where we employ the boundary con- 
dition Ai{x^y) = at the circumference R. 
The superfluid current density is given by ji = 
^^Im [A*djAj + A*d,Aj + A*djAi] . We use known 
GL parameters appropriate for our settings (lol. [HI . [l^ . 
The employed GL parameters [is} are for P = 3.2 MPa, 
with a sufficiently strong H; H » Hd {H = 26.7 mT 
in Q). These are precisely known and well documented 
thanks to the long history of this research field[l,[l]- The 
external rotation vector f2 is parallel to the cylindrical 
axis z. For t = T/T^ > 0.75 (T^ = 2.4 mK) the A-phase 
is stable. 

We now consider the stable state at rest. Our OP 
is the chiral axial states j4±, which are degenerate at 
f2 = 0. The external rotation removes this degener- 
acy. We denote a state's winding number combination 
as {w+,w-), where A± = |A±|e™±^ [9 is the azimuthal 
angle) . Under our convention for rotation sense (counter 
clockwise rotation) the (0,2) state is stabilized over the 
(2,0), that is, the dominant component is A^ with the 
winding w+ = 0, as shown in Fig. 1. A^ with the 
winding W- = 2 is induced at the boundary because the 
gradient terms such as A+V'^A- inevitably induce other 
components at places where the main OP is spatially var- 
ied, at the boundary and vortex core in our case. Under 
counter clockwise rotation the mass current flows less for 



(0,2) than for (2,0). Thus at rest and in lower rotation 
regions, (0,2) is the ground state. In Fig. 1 we display 
the OP profiles. The main component A+ has no phase 
winding, while the minor component j4_ has 47r phase 
winding, which leads to spontaneous mass flow around 
the outer boundary even at rest. 

We now examine a stable vortex under rotation. There 
are a variety of possible vortex states. The general form 
of the winding number combination in the axis symmetry 
confined in a cylinder is given by {w+ = n, w- = n ± 2) 
with n an integer, meaning that through the gradient 
term, two components with w+ = n and w_ = n ± 2 are 
coupled. This includes (-1,1), (1,-1), (0,2), (2,0), (1,3), (- 
3,-1), (2,4), etc. Note that in general, vortices with higher 
winding number are less plausible. The (-1,1) state is 
comparable in energy with the (1,3) vortex. The former 
is more stable than the latter for systems with smaller 
radii, R — 2.5, 5.0, 10 fim in the lower rotation region. 
However, since the (2,0) state in that region has the low- 
est energy, (-1,1) is not realized. The (1,-1) vortex is 
another candidate, but it is not stable in our calculation. 
The higher winding vortex states, such as (2,4) become 
more stable than the (1,3) vortex at higher rotations. 
From these arguments, which we have confirmed by de- 
tailed numerical calculations, we are left with only one 
stable state (1,3) next to (0,2) upon increasing rotation. 
Figure 2 displays the spatial structures for two compo- 
nents and phase windings of the (1,3) vortex where the 
dominant and minor components are A^ and A- with 
w-j- — 1 and = 3, respectively. 

For (0,2) versus (1,3), it is obvious that under rotation 
the energy gain for (0,2) due to j • 17 is less than that for 
(1,3), which leads eventually to a phase transition from 
(0,2) at rest to a single vortex with (1,3) with increas- 
ing fi. We have performed extensive calculations to find 
the critical rotation speed f2c as a function of R that 
allows us to safely locate the parameter regions in the 
realistic case of R = 1.5 mm. A difficulty of the numer- 
ical calculations is that we must treat two quite distinct 
length scales, the coherence length ^ '--^ 10 nm and the 
system size i?=1.5 mm, on an equal footing in the nu- 
merics. The possibility of a singular core of ^ scale was 
neglected in previous calculations. Here we have solved 
the full GL functional to take the singular core into ac- 
count, using fine mesh points in the numerics. Buchholtz 
and Fetter [ij] pointed out the possibility of a singular 
vortex in a narrow channel. 

We show the results of our calculations in Fig. 3, show- 
ing the size R dependence of the critical rotation speed 
flc- We have computed various sizes up to i? = 100 /um at 
a fixed temperature t = T/Tc — 0.95. The inset shows an 
example determination of flc for R = 0.8 ^m. It can be 
seen in the main panel in Fig. 3 that flc = 0.0 6rad/sec 
for t~ 0.95 at i? = 1.5 mm, estimated by an extrapo- 
lation. It should be feasible to realize this relationship 
experimentally. 
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FIG. 1: (Color online) Spatial structures of the order param- 
eters \A+\ (left hand side) and \A-\ (right hand side) for the 
(0,2) state. The sense of the phase winding and its number 
are shown. R — 0.8 /im and Q — 100 rad/sec. 




FIG. 2: (Color online) Spatial structures of the order param- 
eters, l^+l (left hand side) and \ A-\ (right hand side) for the 
(1,3) state. R = 0.8 fim and Q — 100 rad/sec. 



For the temperature dependence of f2c we have checked 
that the critical speed flc slightly increases with lowing 
T. We also found that flc delicately depends on R. By 
using an appropriate R we can tune flc to bring it to an 
experimentally available rotation region. A maximum 
rotation speed of fl = 12 rad/sec can be achieved at 
ISSP, University of Tokyo. Thus, according to Fig. 3, R 
must be above i? > 100 /im. 

Having found that the singular vortex characterized by 
(1,3) is stabilized under rotation above flc, we analyze 
the spectral structure of the core localized bound states 
in the (1,3) vortex with the 3D Fermi surface. The sys- 
tem is three dimensional. We solve the following BdG 
equation under these conditions to examine the Majo- 
rana zero energy state: 



H{r) 
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ff(r)--^ 
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'±j:±{A±{r)n± + in±A±(r)} and 
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-Ep-VtL^. We set 



:^e±^e(#±^- 



and Lz = —ih-^ in the cylindrical coordinates r =^ 
{r,6,z). Equation ^ is self-consistently coupled with 
the gap equation A±{r) ^ ^T,\E^\<Ej'"q'^'^)OTM^) " 
Uq{r)a^v*{r)]fq, where = l/{e^i/''BT + i) the 
Fermi distribution function at a temperature T and we 
henceforth set r = 0. The cutoff energy in the gap equa- 
tion is set as Ec = Ep and the coupling constant g is cho- 
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FIG. 3: (Color online) Size 7? dependence of the critical angu- 
lar velocity fie from (0,2) at rest to (1,3) at a higher rotation 
(T/Tc = 0.95). An extrapolated value of f2c = 0.06 rad/sec is 
found at R= 1.5 mm. The inset shows the energy differences 
(arbitrary scale) between the two states as a function of Q. for 
R = 100 tiTCi. 



sen to be A:fC = 2(A/£:f)~^ = 14.3 with A = max |^±(r)|. 

We assume cylindrical symmetry of the OP having a 
single vortex at r = 0, i.e., A±{r) — A±{r)e™^^ with 
winding numbers wj^ = n and W- = n + 2. Also, we 
impose a periodic boundary condition with periodicity 
Z — 250k along the z axis and a rigid boundary con- 
dition at r = R = 250kp^ in Eq. ([4]), which allows us 
to classify the eigenstates of Eq. (jl]) with the azimuthal 
quantum number G Z and the linear quantum num- 



be 



— with 



e Z. To this end. 



we find 
It 



Uq(r)=Uq{r)e"io'^e"i'\ Wq(r) = Uq(r)e'(9«-"- 
is worth mentioning that for the (1,3) vortex the eigenen- 
ergy under arbitrary il yields a linear shift from the value 
without rotation Eq{^ — Q) as 



Eq=Eq{n = {))- 



1 



(5) 



under low rotation within ^l-^Ep. 

A singular vortex produces a quantum well potential 
for quasiparticles around the core within r <g; ^, which 
gives rise to a low energy eigenstate bounded at the core 
via Andreev scattering. Under no rotation, the eigenen- 
ergy of the core bound state in the 3D singular vortex 
with arbitrary winding n can be analytically obtained 
from Eq. ^ within Igel^^F^ and \qz\<.kF 
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1 



sin(Q!)a;i, (6) 



where sin^(Q!) = 1 



and m G Z. Here, 



assuming 



A+(r)=: Atanh(r/^) and A_(r) = 0, we find ujo~—n 
and ~ A. The remarkable feature of Eq. ([6]) is that 
the eigenenergy can become exactly zero in a vortex of 
odd number n, in spite of the discreteness of the energy 
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due to the narrow core. Also, the zero energy state is 
dispersionless on the hnear momentum . Note that the 
zero energy state appears at = 1 in the (1,3) vortex 
with n— \ and at = in the (—1,1) vortex with n = — \. 
It is found from Eq. ([5]) that the eigenenergy always lies at 
zero energy against the increase of the rotation frequency 

n. 

To clarify the stability of the Majorana zero energy 
state, we present in Fig. [4] the low-lying quasiparticle 
states of the (1,3) vortex as functions of qe and qz, ob- 
tained from a hill numerical calculation of Eq. (jj]) along 
with the gap equation. Note that |^+(r)| ^ |A_(r)|. It 
can be seen that the zero energy state appears at qg — l, 
which is dispersionless situated exactly at i? = within 
\(lz\ < kp. In Fig. m^b), the fact that the bulk excita- 
tion gap changes from the maximum value \Eq\ = A at 
qz — to \Eg\ a.t qz — zkkp indicates the presence of 
point nodes, due to the A phase, situated at the north 
and south poles on the 3D Fermi surface. Since the zero 
energy state results from the resonance between two low- 
lying modes, such as the core bound and edge modes, it 
yields two-fold degeneracy and its wave function is local- 
ized at both the vortex core and edge, where Uq{r)^kpr 
inside the core region r <C ^. It should be emphasized 
again that the flat band structure of the zero energy 
states along qz is independent of the rotation frequency 
fi, because the angular momentum of the quasiparticles 
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FIG. 4: (Color online) Energy spectrum Eq with = Q (a) 
and qe = 1 (b) in the (1, 3) vortex under no rotation f2 = 0, 
where A = max|A±(r)|. The branches labeled "Core" and 
"Edge" in (a) are eigenstates having wave functions tightly 
bounded at the vortex core and the edge, respectively. The 
zero energy state (ZES) appears within \qz \ Kkp at qe = l- 



qe = 1 is canceled by the vorticity +1 and the chirality 
-f 1 of the Cooper pair, as seen in Eq. 

The Majorana core bound state and the edge mode in- 
evitably hybridize, causing level splitting when the sys- 
tem is small enough. For the detection of the Majo- 
rana state, an experiment could be devised to search for 
the energy splitting, which is absent in ordinary vortex 
bound state. The energy splitting could be detected by 
sound attenuation or electromagnetic absorption experi- 
ments. 

It is possible to study similar Majorana physics for a 
different configuration with the field orientation perpen- 
dicular to the plates, as is realized in Ref. 0]. Since the 
d-vector is now free to rotate in the plate plane, the su- 
perfluid becomes "spinful" . By applying a stronger field, 
we can explore the spin polarized Ai phase with a sub- 
stantial temperature window which could be wider than 
the A phase window at iJ = under if ~ 10 t[3|. 

In conclusion, we have proposed an elegant method to 
create the Majorana zero energy state by controlling the 
Z-vector and d-vector of the '^He-A phase via the bound- 
ary condition and external field. The critical rotation ve- 
locity is accurately and realistically estimated as a func- 
tion of the system radius R based on the established GL 
functional form. The spectral properties of the stable 
singular vortex are examined by solving the BdG equa- 
tion for a 3D Fermi surface, explicitly demonstrating that 
there exists a dispersionless Majorana state localized in 
a vortex core. The utility of this 3D Majorana state for 
other fields, such as quantum computation, etc., is open 
for future study. 

The authors thank O. Ishikawa and T. Ohmi for pro- 
viding expertise on superfiuid ''He. 
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